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ABSTRACT 



Growing theoretical evidence suggests that the first generation of stars may have been quite 
massive (~ 100 — 300 Mq). If they retain their high mass until death, such stars will, after 
about 3Myr, make pair-instability supernovae. Theoretical models for these explosions have been 
^ ' studied in the literature for about four decades, but very few of these studies have included the 

, effects of rotation and none ever employed a realistic model for neutrino trapping and transport. 

■ Both turn out to be very important, especially for those stars whose cores collapse into black 
holes (helium cores above about 140 M0). We consider the complete evolution of two zero- 

■ metallicity stars of 250 and 300 M0. Despite their large stellar masses, we argue that the low- 
, metallicities of these stars result in negligible mass-loss. Evolving the stars with no mass-loss and 
' including angular momentum transport and rotationally induced mixing, these two stars produce 

fH ^ helium cores of 130 and 180 Mq respectively. Products of central helium burning (e.g. primary 

nitrogen) are mixed into the hydrogen envelope, which can dramatically change the expansion of 
the envelope, especially in the case of the 300 Mq model. Explosive oxygen and silicon burning 
^ • cause the 130 Mg helium core (250 Mq star) to explode, but explosive burning is unable to drive an 

' explosion in the 180 Mq helium core and it collapses to a black hole. For this star, the calculated 

angular momentum in the presupernova model is sufficient to delay black hole formation and 
the star initially forms a 50 Mq, 1000km core within which neutrinos are trapped. Although 
the star does not become dynamically unstable, the calculated growth time of secular rotational 
instabilities is shorter than the black hole formation time, and such instabilities may develop. 
' The estimated gravitational wave energy and wave amplitude would then be Eqw ~ 10~^ Mq 

and w 10^^^ /d{Gpc), but these estimates are very rough and depend sensitively on the non- 
linear nature of the instabilities. After the black hole forms, accretion continues through a disk. 
The mass of the disk depends on the adopted viscosity, but may be quite large, up to 30 Mq 
when the black hole mass is 140 Mq. The accretion rate through the disk can be as large as 
1-10 MqS~^. Although the disk is far too large and cool to transport energy efficiently to the 
rotational axis by neutrino annihilation, it has ample potential energy to produce a 10^** erg jet 
driven by magnetic fields. The interaction of this jet with surrounding circumstellar gas may 
produce an energetic gamma-ray transient, but given the redshift and time scale, this is probably 
not a model for typical gamma-ray bursts. 

Subject headings: gamma rays: bursts — stars: supernovae, nucleosyntliesis 

1. Introduction stars (Ostriker & Gnedin 1996) contained many 

extremely massive members, from one hundred to 
Simulations of the collapse of primordial molec- several hundred solar masses (e.g. Larson 1999; 

ular clouds suggest that the first generation of Bromm, Coppi, & Larson 1999; Abel, Bryan, 
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& Norman 2000). While details of the mass 
function and the interaction of these stars with 
their environment have yet to be worked out, 
up to 1 % of the baryonic mass of the universe 
might have participated in this generation of stars 
(Abel, private communication). Such stars (M 
k. 100 Mq) will reach carbon ignition with he- 
lium core masses in excess of about 45 Mq and 
will encounter the electron-positron pair instabil- 
ity, igniting carbon and oxygen burning explo- 
sively (e.g., Barkat, Rakavy, & Sack 1967; Woosley 
& Weaver 1982; Bond, Arnett & Carr 1984; Carr, 
Bond, & Arnett 1984; Glatzcl, El Eid, & Fricke 
1985; Woosley 1986; Heger and Woosley 2000). 
If explosive oxygen burning provides enough en- 
ergy, it can reverse the implosion in a giant nu- 
clear powered explosion. As the mass of the he- 
lium core increases, so does the strength of the ex- 
plosion and the mass of ^^Ni synthesized. Indeed 
masses of ^^Ni ejecta over 40 and explosion 
energies approaching 10^^ erg are possible, with 
light curves brighter than 10'*'*crgs~^ for several 
months. Truly these are "hypernovae" (Woosley & 
Weaver 1982; Paczynski 1998). But they have not 
hitherto been associated with gamma-ray tran- 
sients (GRTs) because no relativistic matter is 
ejected. 

However, going to still more massive helium 
cores (over 140 Mq) a new phenomenon occurs as 
a sufficiently large fraction of the center of the star 
becomes so hot that the photodisintegration insta- 
bility is encountered before explosive burning re- 
verses the implosion. This uses up all the energy 
released by previous burning stages, and, instead 
of producing an explosion, accelerates the collapse. 
A massive black hole is born inside the star. This 
has been known for some time, but here we con- 
sider what might occur if, in addition to its high 
mass, the helium core is endowed with a moder- 
ate amount of rotation. We follow zero mctallic- 
ity stars from the main sequence where they are 
assumed to rotate rigidly with a ratio of surface 
centrifugal force to gravity of 20 % - comparable 
to what is seen in 0-stars today. By combining the 
angular momentum transport physics described in 
Hegcr, Langer, & Woosley (2000) up until stellar 
collapse with the two-dimensional rotating core- 
collapse code developed by Fryer & Heger (2000), 
we are able to follow the entire life and death of 
this rotating 300 Mq star. Bond et al. (1984) dis- 



cussed the effects of rotation on the collapse of 
these massive stars and we compare our simula- 
tions with their estimates for neutrino emission 
and gravitational wave emission. 

Because of their high angular momentum, such 
massive helium-depleted cores do not immediately 
collapse to black holes, but instead form a hot 
dense, neutronized core of '^SOMq which contin- 
ues to accrete matter for nearly 1 s before collaps- 
ing inside its event horizon. Although this "proto- 
black hole" is not dynamically unstable, it rotates 
sufficiently rapidly to develop secular instabilities 
that may cause it to break into multiple cores, each 
of which collapses to a black hole and then merges 
with the rest. Even without this fragmentation, 
the collapse of a 300 Mq star turns out to be a 
strong source of gravitational waves with energies 
and wave amplitudes of: Eqw ~ 10~^ Mq and 
h+ » lO-'^y d{Gpc). 

After black hole formation and following the 
collapse of roughly another 100 Mq, the outer lay- 
ers of the helium core have sufficient angular mo- 
mentum to hang up in a disk which accretes onto 
the ~ 140 Mq black hole at rates in excess of 
1MqS~^. This black hole accretion disk system 
may produce a new class of collapsar-like (Mac- 
Fadyen & Woosley 1999; MacFadyen, Woosley, & 
Heger 2000) GRTs. Although neutrino annihila- 
tion is too inefficient to produce polar outflows 
in such large stars, if magnetic fields can extract 
^ 1 % of the rest mass of the accreted material 
and focus it into ^ 1 % of the sky, the equivalent 
isotropic energy of these GRTs could approach 
10^*5 erg. 

2. Pre-Collapse Evolution of 250 and 
300 Mq Stars 

2.1. Population III Stars 

It has long been thought that the first genera- 
tion of stars after the Big Bang might be charac- 
terized by an initial mass function skewed to more 
massive stars (e.g.. Silk 1983; Carr & Rees 1984). 
More recently Abel, Bryan, & Norman (2000) have 
studied the formation and fragmentation of pri- 
mordial, zero metallicity molecular clouds using a 
three-dimensional code with adaptive mesh. Dark 
matter dynamics, hydrodynamics, and the rele- 
vant chemical and radiative processes were fol- 
lowed down to a scale of 0.5 pc. They concluded 
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that the typical mass of the first-generation stars 
is ~ 100 Mq. Based upon the amount of energy re- 
leased by hydrogen fusion, the need to re-ionize a 
significant fraction of the universe prior to red shift 
5 using the light of such stars, and a crude esti- 
mate of the ionization efficiency of a given ultravi- 
olet photon, one estimates that roughly 0.01 % to 
1 % of the baryonic mass of the universe was incor- 
porated into such stars. For a 10^^ M0 galaxy this 
corresponds to ^ 10"* to 10^ massive stars, an es- 
timate that probably exceeds the rates of merging 
double neutron star binaries or black hole binaries 
in the Milky Way, 10"* - 10^ per lOGyr (Fryer, 
Woosley, & Hartmann 1999). There is thus am- 
ple theoretical basis for assuming the formation of 
an appreciable number of primordial stars having 
mass up to several hundred Mq. 

2.2. Mass Loss 

Because of their short lifetime, there arc no 
nearby examples of these zero-metallicity, high- 
mass stars formed in the first generation of stars. 
Indeed, because the initial mass function of pop- 
ulation I and II stars is skewed much more to- 
ward low-mass stars, there are few examples of any 
nearby massive stars. Unfortunately, this makes 
determining the mass loss of these stars from ob- 
servations very difficult. 

Mass loss in such heavy stars today is domi- 
nated by radiative processes (Appenzeller 1986) 
and the stars probably lose most of their mass 
before dying. For zero mctallicity, however, as 
would have characterized these first stars, radia- 
tive mass loss is probably negligible (Kudritzki 
1999). However, one must still consider mass loss 
driven by nuclear pulsations. The critical mass 
for the onset of the nuclear pulsational instability 
(Schwarzschild & Harm 1959) is very uncertain 
(90 . . . 420 Mq for stars of solar metallicity; Ap- 
penzeller 1986), and the resulting mass loss rate, 
even more so. For very low metallicity, both the 
temperature dependence of hydrogen burning and 
the gravitational potential at the surface of the 
star would be different. AppenzcUer's (1970) es- 
timate that a 130 Mq star would lose ~ 30 Mq 
during its main sequence lifetime is probably an 
overestimate for Pop HI stars. What the value 
would be for stars of 200 and 300 M© is unknown, 
but preliminary calculations by Baraffe, Heger, & 
Woosley (2000) suggest that it could be small. 



These same new calculations also suggest that, at 
solar metallicity, opacity and ionization instabili- 
ties are more important in driving mass loss than 
hitherto realized. 

2.3. Presupernova Evolution of 250 and 
300 Mq, Z = stars 

We thus consider the evolution of two massive, 
rotating. Population III stars of 250 and 300 M© 
evolved at constant mass. By "Population III", 
we mean that the initial composition is 76 % H 
and 24% He with no initial abundance of any- 
thing heavier. The evolution of these stars is fol- 
lowed using a one-dimensional hydrodynamic stel- 
lar evolution code, KEPLER (Weaver et al. 1978), 
from central hydrogen ignition until core collapse. 
Stellar rotation is included as described in Heger 
et al. (2000). That is, mixing and redistribution of 
angular momentum induced by convection, shear, 
Eddington- Sweet circulation, etc., arc followed, 
but the centrifugal force terms are not included 
in the stellar structure equations. Semiconvection 
is treated as in Woosley & Weaver (1995), but no 
overshoot mixing is assumed. So long as the ra- 
tio of centrifugal force to gravity remains small, 
as it did during the epoch followed by KEPLER, 
this should be an accurate approximation. For the 
initial angular momentum of the rotating ZAMS 
stars we assumed rigid rotation with a velocity 
20 % Keplerian at the equator, comparable to the 
typical value of observed 0-stars. 

More details of the evolution of these stars and 
many other Population III models will be given 
elsewhere (Heger & Woosley 2000). Here we note 
only those characteristics relevant to forming black 
holes and producing jets and gravity waves. Af- 
ter central helium depletion, (T = 5xlO^K) the 
masses of the helium cores in the 250 and 300 Mq 
models are 130 and 180 M©, respectively. The hy- 
drogen envelopes are also appreciably enriched in 
helium and products of helium burning with X = 
16.0%, Y = 68.0%, Z = 16.0% for the 250 Mq 
model and X = 19.4%, Y = 72.2%, Z = 8.4% 
(mostly CNO) for the 300 Mq model. The lumi- 
nosities of the two stars are 3.8xl0^''ergs~^ and 
5xl0*°ergs~^, respectively, and the central car- 
bon abundance, 6.0% and 4.8%, by mass. 

An interesting characteristic of these two stars 
is their production of primary nitrogen. As a con- 
sequence of both the initial CNO-deficiency and 
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helium-rich nature of their envelopes (resulting 
from rotationally induced mixing), the hydrogen 
burning shell is characterized by a relatively shal- 
low entropy gradient. Thus the envelope remains 
compact and rapidly rotating with only a small en- 
tropy barrier separating it from the helium core. 
During the growth of the helium-burning convcc- 
tive core, traces of carbon and oxygen are mixed 
into the hydrogen-burning shell by meridional cir- 
culation. An increase of CNO mass fraction to 
as little as > 10~* is sufficient to significantly 
increase the nuclear energy generation rate and 
make the hydrogen burning shell convectivc. More 
significantly, during late stages of helium burn- 
ing the shear between the hydrogen shell and the 
core becomes large enough to lead to a significant 
dredge- up of the helium core, i.e., helium and large 
amounts of helium burning products - carbon, oxy- 
gen, and neon - are mixed into the envelope. Be- 
cause of the short time remaining in the life of the 
star at this point, only part of the carbon and oxy- 
gen arc processed into nitrogen. For the 250 M0 
star, this major mixing event happens at a central 
helium mass fraction of ~25 % while in the 300 Mq 
star it occurs only at the end of central helium 
burning. In the 250 star, the mass fraction of 
primary nitrogen in the envelope is 7.75% (i.e., a 
total of 9.48 M0), while carbon and oxygen have 
mass fractions of 0.26% and 7.97%, respectively. 
Since in the 300 Mq star the dredge-up occurs only 
towards the end of central helium burning, the ma- 
terial mixed into the envelope shows a clear signa- 
ture of elements produced at the end of hot helium 
burning. For the 300 star, the mass fractions 
are: nitrogen - 1.28% (i.e, a total of 1.56 M©), 
carbon - 0.033%, oxygen - 5.68%, neon - 0.964%, 
and magnesium - 0.395%. 

The increased energy generation afforded by 
the (suddenly) larger abundance of CNO leads to 
the formation of red supergiants in both stars. 
The radius of the 250 model expands from 
2.3 X 10^^ cm to ~10^^ cm by the time it encounters 
the electron-positron pair instability; the 300 M0 
model, from 4x10^^ cm to ~1. 5x10^** cm. The 
mixing events occur ~10^ yr before core collapse in 
the 250 M0 model and ^10'* yr before collapse in 
the 300 Mq model. Note that all these estimates of 
mixing, nitrogen nucleosynthesis, and time-scales 
are sensitive to the highly uncertain physics of con- 
vection and rotational instabilities. 



The new large radius and high metal content in 
this supergiant phase can cause significant mass 
loss and may remove the envelope, particularly if 
the star is in a binary. Once the helium core is 
revealed, unless a small amount of hydrogen re- 
mains, an even more pulsationally unstable situ- 
ation may be created. Appenzeller (1986) sum- 
marizes evidence that bare helium cores above 
~16M0 may be unstable to nuclear driven pulsa- 
tions. However, especially for the 300 model, 
the time remaining until the star dies is very short. 
If the core can avoid losing an additional 40 Mq 
of hehum (i.e., M & O.OO4M0 y~^), a black hole 
will form. 

3. Collapse and Black Hole Formation 

Following helium depletion, both the 250 and 
300 Mq stars encounter the electron-positron pair 
instability. The 250 Mq star experiences a very 
deep bounce, penetrating into the pair-unstable 
region, then violently rebounds because of the 
excess energy created by explosive oxygen and 
silicon burning. The peak temperature and 
density during the bounce are 6.37 x 10^ K and 
1.08x10'^ gcm~^ (Fig. 1). So much of the core is 
heated to temperatures in excess of ~ 5xlO^K 
that 43.0 Mq of ^^^Ni is synthesized. We esti- 
mate that this is very nearly the maximum mass 
that can explode by nuclear burning alone and 
so represents a nearly maximal mass of ^^Ni and 
explosion energy, 9 x 10^^ erg of kinetic energy at 
infinity, almost 100 times that of an ordinary su- 
pernova. Heger et al. (2000) have calculated 
a peak luminosity for this model of ~10^^ergs~^ 
(Mboi = —21) lasting for about 150 days. Roughly 
the same peak luminosity and duration charac- 
terizes the event with and without the hydrogen 
envelope, but in the case of the helium core the 
peak is delayed about 100 days and completely 
powered by ^^Co decay. Since the kinetic energy 
of the explosion is so high, the interaction with 
circumstellar matter might also be quite brilliant. 
These are truly "hypernovae" (Woosley & Weaver 
1982), but no strongly relativistic matter is ejected 
and this explosion will not produce a GRT. 

The 300 Mq model (180 M0 hehum core) on the 
other hand is so tightly bound and gains so much 
kinetic energy prior to oxygen burning that even 
the fusion of the entire core to silicon and iron is 
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unable to reverse its infall. It must make a mas- 
sive black hole. At the latest time reliably cal- 
culated with KEPLER, the central temperature 
had reached 10^^° K and the nuclei in the center 
of the star had photodisintegrated to 72 % alpha 
particles and 28% nucleons {Ye = 0.50). Unlike 
what is seen in ordinary supernova models, even 
though temperature dependent partition functions 
were included in the calculation, the higher en- 
tropy of these models results in complete photo- 
disintegration of heavy elements. The density at 
this point is 5.5xlO''gcm~^, the total kinetic en- 
ergy of infall, 5.8x10^^ erg, and the net binding 
energy, —7.9x10^^ erg. Figure 2 shows the com- 
position, temperature and density structure, and 
distribution of specific angular momentum at this 
time. In subsequent discussions we define this to 
be the t = collapse model. 

3.1. Collapse without Rotation 

Beyond this point the effects of neutrino trap- 
ping are no longer negligible and the Kepler model 
was mapped into a one-dimensional Lagrangian 
core-collapse code (Herant et al. 1994, Fryer et 
al. 1999) which included the effects of general rel- 
ativity and neutrino transport. It is useful to 
compare the ensuing collapse of this 300 Mq star 
to the core-collapse of an ordinary 15 super- 
nova progenitor. Just like ordinary supernovae, as 
the core of the 300 star contracts, the rate of 
electron capture increases and the removal of elec- 
tron degeneracy pressure and cooling via neutrino 
emission help to create a runaway collapse. How- 
ever, there are several important differences be- 
tween the structures of 300 and 15 Mq stars. 
As noted previously, the entropy in the more mas- 
sive core is larger and thus favors the more com- 
plete photodisintegration of heavy elements and 
alpha particles. General relativity also plays a 
more significant role. The collapse of a 15 Mq core 
halts when its central density exceeds a few times 
lO^^gcm"^, i.e., when nuclear forces and neutron 
degeneracy pressure become important. However, 
the core of a 300 Mq star is so large that it col- 
lapses into a black hole before nuclear forces can 
affect the collapse, less than 1 s after the code link 
at Tc = IQi" K. 

Figure 3 shows the velocity, electron fraction 
(Ye) 5 temperature and density profiles of the 
300 Mq core at times zero, 500, 800, and 930 ms 



after the code link. Before the initial collapse, 
the entropy of the core is ~ 10A;b per nucleon 
(compare to ~ 1 — 2fcB per nucleon for most 
corc-collapsc supernovae). Hence, for a given cen- 
tral collapse density, the 300 Mq star is hotter 
and neutrino emission is more efficient. Neutri- 
nos from electron capture stream out from the 
core, quickly lowering its lepton number and low- 
ering the electron degeneracy pressure. However, 
because the core of a 300 Mq star is so much big- 
ger than a 15 Mq supernova progenitor, electron 
neutrinos are trapped at a lower central density 
(Bond et al. 1984). Initially, neutrino absorp- 
tion dominates the electron neutrino opacity, but 
as the temperature rises, pair scattering quickly 
becomes the most important opacity source. Al- 
though /X and r neutrinos escape the collapsing 
star more easily, trapping the electron neutrinos 
halts the dclcptonization in the core. But at this 
point, the collapse of the core is already inevitable. 
The large population of electron neutrinos in the 
core actually causes the central electron fraction 
to rise, as the core strives to reach an equilib- 
rium between the electron/anti-electron neutrinos 
and the electron/proton fractions in the core (see 
Fig. 3). 

In our simulations, we follow the collapse until 
a part of the star falls below its last stable orbit 
(for the non-rotating star this is 6GMenciosed/c^)- 
We use the simplifying assumption that the core 
"forms a black hole" at this point. A rigorous de- 
termination of the first trapped surface would re- 
quire a multi-dimensional general relativistic cal- 
culation that is beyond the scope of this paper. 
The exact point in the core where the initial black 
hole forms is also very sensitive to the neutrino 
transport (especially upon the /U and r neutrinos 
which dominate the neutrino cooling) . In our sim- 
ulations, the black hole first forms at ~ 20 Mq, 
but the entire inner 35 Mq of the core is very 
close to its last stable orbit and relatively small 
changes in the n and r transport could give an 
initial black hole mass anywhere in the range of 
15-35 Mq. However, no matter where the black 
hole initially forms, without rotation, most of the 
helium core will very quickly become part of the 
black hole. 

The bulk of the potential energy released dur- 
ing the collapse is dragged into the black hole. 
Without rotation to drive a dynamo, it is unlikely 
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that a strong magnetic field will emerge during 
the collapse and magnetic fields will not be able 
to tap any of the potential energy. Likewise, a 
non-rotating collapse will not emit much energy 
in gravitational waves. Bond et al. (1984) ar- 
gued that much of the potential energy could be 
released in a neutrino fireball. However, their an- 
alytic calculation overestimated the temperature 
at the neutrinosphere, which drastically overesti- 
mated the total neutrino luminosity. In our sim- 
ulations, most of the neutrinos are trapped in the 
flow and only 1% of the gravitational potential 
energy is released in neutrinos, most in /i and t 
neutrinos (Fig. 4). The collapse of a non-rotating 
300 Mq star proceeds without so much as a whim- 
per. 

3.2. Collapse With Rotation 

However, these stars are rotating and this plays 
an important role in the collapse. As with the 
non-rotating model, the rotating collapse is simu- 
lated by first mapping the KEPLER output into 
the one-dimensional neutrino transport code. The 
effect of rotation is included by adding a cen- 
trifugal term, acent = j'^/i"^, to the force equa- 
tion and conserving angular momentum locally, 
j = j{m,t) = j{m), for each mass m. We map 
j from the KEPLER output (Fig. 2) by setting 
j{r) = rri(r) and add the equatorial force to the 
entire star. This approximation obviously over- 
estimates the eff'ect of angular momentum every- 
where except in the equatorial plane, but allows us 
to follow the first 1.5 s of the star's collapse quickly 
and gives approximately the correct composition 
and mass distribution. Recall that in the non- 
rotating simulation, a black hole formed 1 s after 
mapping from the KEPLER output into our ID 
neutrino transport code. In the rotating simula- 
tion, centrifugal forces become extremely impor- 
tant (Fig. 5) and slow the collapse. After 1.5 s, 
the central density of the rotating core is only 
5xl0^°gcm~^. 

When the centrifugal force exceeded 25% of 
gravity for any lagrangian zone (Fig. 5), the one- 
dimensional collapse calculation was halted again 
and remapped into a two-dimensional code based 
upon the smooth particle hydrodynamics prescrip- 
tion (Herant et al. 1994, Fryer & Heger 2000). 
Use of the 2D code was delayed until this point 
in order to make the problem computationally 



tractable. The rotational axis was aligned with 
the code's axis of symmetry and the angular mo- 
mentum of each particle was given by the rota- 
tion rate of the zones in the one-dimensional cal- 
culation: j(r) = Q{r)r cos{9) where 9 is the an- 
gle above the plane of rotation. The net effect is 
to give (initially) spherical shells constant angular 
velocity and a total angular momentum equal to 
what they had in the Kepler model at t = 0. 

A total of thrc;e two-dimensional simulations 
were performed this way, each employing a dif- 
ferent treatment of the angular momentum. The 
models will be referred to by three names: A) 
no rotation (this is inconsistent with the assumed 
density structure, but calculated just for compari- 
son); B) the star rotates with the angular momen- 
tum distribution at the end of the ID calculation 
- which is the same as prescribed by the KEPLER 
calculation; in the continued evolution the angular 
momentum of each particle is conserved; and C) 
same as B, but with angular momentum transport 
between particles mediated via an alpha-disk pre- 
scription (adisk = 0.1: see Fryer & Heger 2000 for 
details). In the remainder of this section, we dis- 
cuss the collapse of the star and the resultant black 
hole formation under these three assumptions. A 
summary of results is given in Table 1. 

As might be expected from our non-rotating 
collapse simulation, Model A collapses directly 
into a black hole 250 ms after the beginning of the 
two-dimensional simulation. As before, we assume 
that the star has made a black hole when matter 
at any radius lies within its last stable orbit for 
the given distribution of matter. At that point, 
all particles inside this radius are removed and re- 
placed by a perfectly absorbing inner boundary 
condition at that radius. This way, we can fol- 
low the continued accretion into the black hole 
(Sect. 5). 

The models which include the KEPLER angu- 
lar momentum distribution (Fig. 2) evolve very 
differently from the non-rotating case. As the 
star collapses, neutrinos become trapped in the 
flow, maintaining ^ 0.35. The high entropy 
of the 300 Mq star, along with the rapid rotation, 
halts the collapse of the core and even produces 
a weak "bounce" at a central density of only a 
few times lO^^gcm"^ (Fig. 6). The structure of 
the resulting core is similar to the proto-neutron 
star cores encountered in ordinary core-collapse 
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supernovae, but here it is much larger. UsuaUy, 
in lower mass supernovae, after nuclear forces and 
neutron degeneracy pressure halt the collapse, the 
core bounces, sending a shock into the star. When 
the bounce shock stalls, it forms a proto-neutron 
star core capped by the accretion shock produced 
as the remainder of the star falls onto the core. 
In the case of core-collapse supernovae, the proto- 
neutron star has a mass of roughly IM© and a 
radius of ~ 400 km. In the present simulations 
of the 300 Mq star, our weak "bounce" yields a 
50 Mq "proto-black hole" with a radius of 1000 km 
(Fig. 7). By proto-black hole, we mean a hot dense 
neutronized core that has not contracted inside its 
event horizon, but which would do so if lost its in- 
ternal entropy. 

Material continues to accrete onto this core 
through a shock at a rate of 10-100 MqS^^. By 
1.1s after the beginning of the two-dimensional 
simulations (2.6 s after the initial collapse in 
Fig. 2), the proto-black hole mass has increased to 
nearly 70 (Fig. 8). Neutrino emission, domi- 
nated by n and r neutrinos (since the electron neu- 
trinos are trapped), cools the proto-black hole and 
the radius of the accretion shock slowly shrinks. 
From Fig. 9, we see that the total neutrino lumi- 
nosity is well above lO^^ergs"^ during much of 
the simulation. Such high neutrino luminosities 
suggest the possibility of a supernova-like explo- 
sion or perhaps a GRB (Fuller & Shi, 1999). 

However, inside the proto-black hole, scatter- 
ing and absorption dominate the neutrino opac- 
ity and beyond the ~1000km neutrinosphere, the 
density of neutrinos is too low to allow significant 
numbers to annihilate. The neutrino density is 
roughly 10'^'' ergs^^ cm^'', over an order of mag- 
nitude lower than the lowest neutrino density con- 
sidered by Popham, Woosley, & Fryer (1999) for 
black hole accretion disks. In addition, the bulk 
of the neutrinos emitted by the 300 Mq model are 
a) traveling radially and the neutrino cross-section 
is only large for head-on collisions and b) /i and 
T neutrinos which have lower neutrino annihila- 
tion cross-sections. Combining all of these effects, 
neutrino annihilation injects less than 10^^ ergs~^ 
into the infalling material. Given that the binding 
energy of the material piling up at the accretion 
shock is over 10^^ erg, we conclude that neutrino 
annihilation does not drive an explosion. Neutri- 
nos interacting with nucleons and scattering on 



electrons were also included in the calculation and 
these too provided totally insufficient energy to 
reverse the infall. 

However, shortly after the black hole forms, the 
accretion shock moves inwards and angular mo- 
mentum becomes more important. An accretion 
disk will eventually form and that may produce 
an explosion. Thus, it is important to determine 
when the black hole forms. As the proto-black hole 
contracts, its inner core is slowly compressed down 
to the event horizon. As with the non-rotating 
model, we follow the collapse until the material is 
compressed down to the radius of the marginally 
stable circular orbit (e.g., Shapiro & Teukolsky 
1983): 

rxns = Mbh {3 + ^2 - V(3-^i)(3 + ^i + 2Z2)} 

(1) 

where 

(2) 

Z2 ^ ^SaWMin + ^f, (3) 

and flBH/^'^BH and Mbh arc the dimensionless an- 
gular momentum and mass of the black hole. At 
any given radius in the proto-black hole, we calcu- 
late both flBn/AfBH and A^bh and can thus deter- 
mine when that radius falls below the marginally 
stable orbit. When this happens, we stop the sim- 
ulation and remove all of the particles within that 
radius. The total mass of these particles forms 
the initial black hole mass. In this manner, al- 
though we do not physically model the collapse to 
a black hole, we get roughly both the correct time 
of collapse and the correct initial black hole mass 
(Table 1). 

Or do we? Recall that our collapsing star has 
a large amount of angular momentum. As it 
collapses, local conservation of angular momen- 
tum causes it to spin up. As long as the rota- 
tional energy is less than ~ 14% of the gravita- 
tional potential energy, the proto-black hole is sta- 
ble against triaxial deformation. However, when 
T/|W^| > 0.14, instabilities can occur (e.g. Shapiro 
& Teukolsky 1983). By plotting the energy ratios 
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(T/|Vl^|) as a function of radius for the proto-black 
hole 0.5 s, and Os before black hole collapse, we 
see that much of the rotating proto-black hole is 
unstable to triaxial deformations (Fig. 10). The 
growth time for such instabilities is (Schutz 1983): 

T/\W\ f^r\~^ 



where 17 is the rotational velocity at radius r. Fig- 
ure 10 also shows the instability growth time ver- 
sus mass 0.5s and Os before the collapse of the 
proto-black hole. According to our analysis, in- 
stabilities would have time of grow in our rotat- 
ing models, possibly forming smaller clumps which 
would collapse and then merge to form the cen- 
tral black hole (Bond et al. 1984). By com- 
paring the growth time to the time to collapse 
for both our rotating models, we can estimate 
the range of masses in the proto-black hole which 
might develop these instabilities (Table 1). Equa- 
tion (4) is not sufficiently accurate to decisively 
prove that such instabilities will develop before 
the proto-black hole collapses and such verifica- 
tion awaits three-dimensional simulations which 
can model the star over tens to hundreds of or- 
bits. Even if these instabilities occur, they will 
not affect the accretion onto the black hole signifi- 
cantly (Sect. 5) because the individual clumps will 
merge and form a central black hole at roughly the 
same time as our assumed stable model. 

Although the qualitative picture we have de- 
scribed is the same for both rotating models (B 
and C) , the actual formation time of the black hole 
differs considerably. This is because in Model C, 
angular momentum is transported out of the core 
using the a-disk prescription. This lowers the spin 
of the black hole {ubh) which increases the radius 
of the marginally stable orbit (rms) and ultimately 
causes the black hole to form more quickly. By 
choosing a high disk viscosity (adisk = 0.1), our 
Model C estimates the extreme effect of viscosity 
and the true answer probably lies between Mod- 
els B and C. Fortunately, these differences do not 
affect our quantitative estimates of gravitational 
wave emission or GRT energies (Sections 4 and 
5). 



4. Gravitational Waves 

Even if the core does not break into several 
pieces, because r/|T4^| exceeds the secular insta- 
bility criterion (~ 0.14), the star may still de- 
velop a bar-like configuration. Gravitational wave 
emission for rotating bodies in the context of su- 
pernova collapse has been studied in great detail 
(see Zwerger & Miiller 1997 or Rampp, Aliiller, 
& Ruffert 1998 for reviews). We will use here 
the expressions derived by Zwerger & Miiller for 
the quadrupole wave amplitude, the gravitational 
wave field {h+) and the total gravitational en- 
ergy emitted (equations 20-22 of Zwerger & Miiller 
1997). In all of our simulations, l^fg | peaks near 
1.5xl0'^cm, nearly 4 orders of magnitude greater 
than most of the simulations of core-collapse su- 
pernovae by Zwerger & Miiller (1997). The ampli- 
tude of the corresponding waveform is also high: 
h+ = 1.3xlO-27d(Gj3c). Recall that 300 M© 
stars are only likely to form at very high red- 
shifts {z « 15), corresponding to distances beyond 
7.5 Gpc for Ho ~ eOkms'^Mpc'^ For most of 
these stars, h+ « 10~^^. The total energy emitted 
in gravitational waves during the first few seconds 
of collapse is ~ 1O~^M0C^. During this time, a 
100 Mq black hole forms and roughly 10~^ % of its 
rest mass energy is converted into gravitational 
waves. Like the neutrino emission, our calcula- 
tion of the energy lost to gravitational waves is 
much smaller than that predicted by Bond et al. 
(1984). However, if the core of the proto-black 
hole breaks up into pieces, the gravitational wave 
emission could be much larger. 

We have only scratched the surface of our un- 
derstanding of the collapse of these massive stars, 
but already they show some potential as sources 
of powerful gravitational waves. Unfortunately, 
the development of secular instabilities requires 
the ability to model the proto-black hole over 
many orbital periods in three dimensions. While 
we can follow the collapse all the way to black 
hole formation in two dimensions, the develop- 
ment of instabilities (and a more accurate esti- 
mate of gravitational wave emission) awaits fu- 
ture three-dimensional simulations. We would be 
happy to provide any of our models to researchers 
who would like to attempt this problem. 
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5. Black Hole Accretion and Gamma-Ray 

Transients 

Once a black hole has formed, the rest of the 
dense inner core quickly collapses inside, but in the 
rotating models, some material along the equator 
is slowed by centrifugal force. If there is sufficient 
angular momentum to form an accretion disk, this 
material will be slowed and pile up outside the 
hole. The disk that forms would be similar to 
those studied by Popham et al. (1999). Such black 
hole accretion disk systems are believed to power 
gamma-ray bursts, either through neutrino anni- 
hilation or magnetic field powered jets. To assess 
whether the collapse of 300 M© stars might pro- 
duce some sort of gamma-ray transient, we must 
first calculate the disk and black hole properties. 

We follow the accretion into the black hole by 
placing an absorptive inner boundary in the core 
at the marginally stable orbit. As particles move 
inside the marginally stable orbit (rms); their mass 
and angular momentum is added to the black hole. 
The boundary condition corresponding to the last 
stable orbit also expands as the black hole gains 
mass. For the non-rotating model, almost all of 
the potential energy of the accreting material was 
carried into the black hole (recall Fig. 4). Even 
in the rotating models, much of the star collapses 
onto the black hole before a stable accretion disk 
can form. When angular momentum is conserved 
locally (Model B, Fig. 11), a stable disk forms 
only after roughly 140 M© of the star has accreted 
(Fig. 11). 

To produce a GRT, energy generated in the 
disk and transported, either by neutrino emission 
or magnetohydrodynamical processes to matter 
above and below the black hole along the rota- 
tional axis, must generate sufficient pressure and 
momentum to reverse the implosion of material 
accreting along the poles (MacFadycn & Woosley, 
1999). The evolution of the polar accretion and 
disk accretion rates and the efficiency for this 
transport thus dictate if and when an explosion 
is likely to occur. Figure 12 shows the total ac- 
cretion rates for our three simulations along with 
an analytic estimate assuming the mass falls in at 
the free-fall time: 



'0 



2 y GMenclosed 



(5) 



This analytical estimate does not take into ac- 
count pressure forces or angular momentum, both 
of which affect the accretion rate. The accretion 
rate in the non-rotating model shows a large dip 
in the accretion rate over the mass range between 
~ 40 — 80 Mq. This occurs because nuclear burn- 
ing injects energy into the collapsing material and 
slows its infall. This effect is not seen in the rotat- 
ing models because the material in the proto-black 
hole burns into iron elements before the formation 
of the black hole. Later, after black hole forma- 
tion, the material does not have this additional 
reservoir of energy. The accretion rate decreases 
at lower black hole masses in the rotating models 
because some of the material has enough angular 
momentum to support itself in a disk. 

In both of our rotating simulations, a disk forms 
in a wedge extending up to 35-40° above the equa- 
tor. When we include angular momentum trans- 
port (Model C), the angular momentum in the 
outer layers of the helium core increases (as these 
layers gain angular momentum from the core), and 
the disk in Model C is initially larger both in ex- 
tent and mass than in the conserved case (Model 
B). Although the material in the disk falls much 
slower into the black hole than along the poles, the 
density in the disk quickly becomes so much higher 
than in the polar regions that it dominates the ac- 
cretion onto the black hole (beyond a black hole 
mass of 120 Mq, disk accretion makes up > 70% 
of the black hole accretion). By the time the to- 
tal accretion rate drops below 10MqS~^, the ac- 
cretion along the polar region (the entire region 
beyond the 40° disk) accretes less than 3 Mq 
and this rate is dropping rapidly. Unfortunately, 
our limited resolution makes it difficult for us to 
follow the accretion rate when it falls below a few 
M©s~^. In Model B (where the angular momen- 
tum is locally conserved), the accretion rate in the 
disk will drop as abruptly as the rate in the poles 
as the disk becomes entirely supported by centrifu- 
gal forces. However, when we include angular mo- 
mcntimi transport (Model C), the disk will con- 
tinue to accrete at a rate of ~ 1 — IOMqS"^ (for 
an adisk = 0.1) until the entire ~ 40 M© disk has 
accreted onto the black hole. 

It is unlikely that any jet will form before the 
polar regions clear and the accretion rate along the 
poles drops below some critical rate. This critical 
rate depends upon the mechanism used to power 
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the explosion. For the neutrino-driven mechanism, 
the accretion rate in the polar region must drop 
below 1-10 Mqs"^ so that the infalling material 
is optically thin to neutrinos. When the total ac- 
cretion rate falls below IOMqS"^, corresponding 
to an accretion rate in the 50° polar region of 2- 
SMqS"-'^, the conditions around the black hole 
may drive a GRB explosion. At this time, we 
can compare the properties of our black hole ac- 
cretion disk systems with those models produced 
by Popham et al. (1999). The black hole mass 
and spin for both our rotating models are listed in 
Table 1. The total disk mass available for rapid 
accretion is roughly 170 — Mgjj. 

Popham et al. (1999) considered black hole 
masses between 3 and 10 M©. The efBciency at 
which neutrino annihilation converts the gravita- 
tional potential energy released into fireball en- 
ergy drops dramatically as the black hole mass 
increases. As the event horizon increases, a) the 
density and temperature in the disk decreases, de- 
creasing the total neutrino flux and neutrino en- 
ergy and b) the volume of the annihilation region 
increases, reducing the neutrino density. From 
Popham et al. (1999) we find that the neutrino 
annihilation conversion efficiency goes roughly as 
Mg^ where n is roughly 2-3. Even if the disk ac- 
cretes at a rate of 10 Mq s~^, the energy deposited 
by neutrino annihilation will not be high enough 
to drive an explosion. Note that in our simula- 
tions, the accretion rate through the disk in all of 
our simulations is less than IOMqS"^ when the 
polar region has cleared. 

Jets driven by magnetic fields may produce a 
strong explosion though. It is difficult to make 
any quantitative predictions about magnetic-field 
driven GRBs, simply because the exact mecha- 
nism is not well understood. However, by the 
time the total accretion rate has dropped below 
1OM0S^^, a well-defined disk has formed, which 
is required for most of the GRB engines driven by 
magnetic fields. The total potential energy avail- 
able to drive an explosion is egpin emho -^disk , 
where the potential efficiency of the accretion disk 
(fspin) is ^ 0.11 for a/M = 0.74. Assuming 
that the efficiency of our magnetic field driven 
explosion (cmho) is 0.1, the yield from Models 
B and C are 6x10^'^ and 10^^ erg, respectively. 
Beamed into 1 % of the sky, these bursts would 
have inferred isotropic explosion energies of nearly 



Iq56 (^j-g^ easily visible at rcdshift 10 and roughly 
an order of magnitude more energetic than ordi- 
nary collapsars. 

The time scale for these bursts can be crudely 
estimated from the a-disk model (Shakura & Sun- 
yaev 1973), TvIsc ^ /{ctH^flK), with a, the disk 
viscosity parameter, Qk, the Keplerian angular 
velocity, and r and H the radius and thickness of 
the disk, respectively. For r ~ 3H ~ 5000 km 
(Fig. 11) and a ^ 0.1, the time scale is very 
roughly 10 s, but there is considerable uncertainty 
in both a and the other parameters of this equa- 
tion. That this value is approximately equal to the 
duration of many common GRBs is probably co- 
incidental, especially given the large redshift and 
time dilation we expect for 300 Mq stars at death 
(Sect. 2). However, it does suggest the possibility 
of other forms of GRTs with time scales of per- 
haps minutes and a hard x-ray spectrum. Such an 
energetic jet would also likely disrupt the star in 
a gigantic "hypernova" explosion. 

If it does not, then the outer part of the hydro- 
gen envelope will accrete in about a day and has 
enough angular momentum to form a disk (Fig. 2). 
This may make an even longer fainter transient of 
some sort and may be an important source of X- 
ray photons in the early universe. 
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Table 1 
Collapse of 300 Star 



Model'* 


-^Unstable 

(Mo) 


rrt b 

-'Coll 

(s) 


(Mo) 


«bh/^bh 


Afproto-BH 

(Mo) 


^disk 


(M©)^ 


c 


A 




1.75 


13 













B 


5-40 


3.7 


35 


0.73 


90 


3.5 


136 


0.74 


C 


12-35 


2.8 


13 


0.55 


70 


3.6 


122 


0.72 



''Model A neglects rotation. Model B assumes the initial rotation of the progenitor model and conserves 

locally the angular momentum for the; duration of the simulation. Model C uses the same initial rotation as 
Model B, but evolves the angular momcntiun using an a-disk prescription [a^ = 0.1). 

'^Time elapsed after the simulation is mapped from the KEPLER output. 

'^The "disk" quantities (last 3 columns) refer to the time, black hole mass, and black hole rotation at which 
the accretion rate drops below 10 Mq s~^. 
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Fig. 1. — Internal structure of the 250 M© star at maximum central density as a function of mass coordinate. 
Panel a gives the mass fractions of the dominant chemical species. Note that 'iron' denotes the sum of all 
isotopes of the iron group elements. Panel b gives temperature, density, and velocity. Inward movement 
(winfaii) is drawn as dashed line and outward movement (^^expand) as dotted line. The (logarithmic) density 
scale is 3.25 times that of the temperature scale. 
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Fig. 2. — Internal structure of the 300 M© star at a central density of 5.5xl0^gcm~^ as a function of mass 
coordinate. Panel a gives the mass fractions of the dominant chemical species (see Fig. 1). Panel h shows 
temperature, density, and infall velocity (iiinfaii; dashed) at this stage. Panel c gives angular velocity (w) and 
specific angular momentum in the equatorial plane (ieq)- Dashed and dash-dotted lines give, respectively, 




Fig. 3. — Density versus radius of a non-rotating helium core 0, 500, 800 and 930 ms after reaching a central 
temperature of 10^° K and the transition to the general rclativistic code. Rapid neutrino cooling in the core 
causes the inner 20-25 to collapse quickly, decoupling from the rest of the star. It is this inner core that 
first collapses into a black hole. Explosive nuclear burning slows the collapse beyond 40 M©, causing the dip 
in density there. 
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Fig. 4. — Total energy emitted in neutrinos {dotted line: electron, dashed line: /j. and r, solid line: total) 
as a function of time after reaching a central temperature of Tc = 10^" K. As the black hole grows, the 
temperature at the last stable orbit decreases, causing the neutrino emission to decrease dramatically. Note 
that only 1 % of the potential energy released in the collapse actually makes it out of the black hole. 



16 




50 100 150 
Enclosed Mass (Mq) 

Fig. 5. — Ratio of centrifugal force over gravitational force at the end of the KEPLER run (dotted line) 
and 1.5 s later, at the end of the ID collapse simulation. The centrifugal force slows the collapse and the 
density at the end of the ID rotating simulation is only 5x10^° gcm~^. Without rotation, the core would 
have already collapsed to form a black hole. 
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Fig. 6. — Velocity versus radius 250 ms after mapping the simulation into the two-dimensional code (Model 
B). Points correspond to particles at different radii and latitude. Thermal pressure, aided by support from 
angular momentum, slows the collapse of the core, causing a weak "bounce" . The material along the equator 
is supported by angular momentum and collapses much slower than the material along the poles. Hence, 
at any given radius, there is a range of velocities (the slowest speeds correspond to the material along the 
equator). Note the accretion shock developing at 2000km. 
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Fig. 7. — The proto-black hole 0.5 s before black hole formation. Color denotes temperature in 10^ K and 
the vectors represent the direction and magnitude of the particle velocity. The simulation was actually only 
half of this circle and was reflected about the z-axis for display purposes. It used a total of 25,000 particles 
though only ^ 15, 000 are shown in the figure. At this time, the proto-black hole has a mass of roughly 
78 Mq and size of ^ 1100 km. 
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Fig. 8. — Density versus enclosed mass for the collapsing core 1.75 {red; lower curve) and 2.6s {blue; upper 
curve) after reaching a central temperature of Tc = 10^" K. The points correspond to particles at different 
radii and latitudinal angles. The "enclosed mass" is defined by the mass inside a sphere with same radius as 
the particle distance from the center. Because of the range in infall velocities, at any given mass coordinate, 
there is a range of densities. The proto-black hole mass increases as material accretes through the accretion 
shock. However, neutrino cooling causes the core to contract and the accretion shock actually moves inward. 
These two effects cause the density to increase dramatically. 
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Fig. 9. — Neutrino Luminosity as a function of time from Model B. The fi and r neutrinos (dotted line) 
dominate the neutrino emission until black hole formation. Shortly after the black hole forms, the event 
horizon grows beyond the fx and r neutrinosphere (at 2.5 s) and drastically diminishes the neutrino luminosity. 
The electron neutrinos (solid line) do not decrease significantly until the black hole expands enough to 
produce a cool accretion disk. 
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Fig. 10. — T'/IW^I and instability growth time vs. mass in the proto-black hole at two different times for 
Model B: solid line - at black hole formation, dotted line - 0.5s before black hole formation. If T/|Ty| 
is greater than ~ 0.14, the system may develop secular instabilities. From 5-40 M©, the growth time of 
the instabilities is shorter than the actual simulation time, and this region may have time to clump before 
collapse, forming several smaller black holes before merging into a central, spinning black hole. 
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Fig. 11. — The accretion disk 6.5 s after black hole formation assuming local angular momentum conservation 
(Model B). Velocities are color coded in units of 1000 km s^"'^. At this time, the accretion rate through the 
pole is less than 0.1 M0 s~^. At this time, only 5200 particles remain in the simulation (4600 are shown in 
the figure). The color represents radial velocity, showing that a large part of the disk is now stable. 
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Fig. 12. — Total accretion rates for Model A (dotted line), Model B {dashed line), Model C {dot-dashed 
line), and an analytic estimate assuming matter collapses on a free-fall time {solid line). The large dip in 
the accretion rate between 40 and 80 M© in the non-rotating model (Model A) occurs because of nuclear 
burning. Angular momentum transport slows the collapse of material just above the disk in Model C, causing 
the accretion rate to drop slightly quicker in this simulation than in the simulation where angular momentum 
is conserved locally (Model B). Note that the accretion rate drops precipitously once it is below 10 M© s~^ 
and changing the critical accretion rate from 10 to 1 M© s~^ does not change the black hole mass by more 
than 10 %. 
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